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Abstract 

We investigate the magnetic properties of a range of low-dimensional ferromagnets using a com- 
bination of first-principles calculations and atomistic spin dynamics simulations. This approach 
allows us to evaluate the ground state and finite temperature properties of experimentally well 
characterized systems such as Co/Cu(lll), Co/Cu(001), Fe/Cu(001) and Fe/W(110), for different 
thicknesses of the magnetic layer. We compare our calculated spin wave spectra with experimental 
data available in the literature, and find a good quantitative agreement. We also predict magnon 
spectra for systems for which no experimental data exist at the moment, and estimate the role of 
temperature effects. 
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In recent years, spin-polarized electron energy loss spectroscopy (SPEELS) has been de- 
veloped to the point of becoming a powerful method for investigating spin waves at surfaces 
and in thin films. 1 Landmark experiments have included the measurement of the magnon 
spectrum of ultrathin Co films on Cu(OOl), 2 and of a single monolayer (ML) of Fe on 
W(110). ! More recently, an asymmetry in the magnon spectrum of 2 ML Fe on W(110) 
has also been reported - this is a direct signature of the Dzyaloshinskii-Moriya interaction 
present in this system." The experimental accessibility of these properties provides an oppor- 
tunity for comparing the merits of different theoretical approaches used to describe thin-film 
magnets. '~ s Furthermore, the experiments challenge theoreticians to address issues such as 
accurately characterizing exotic ground states arising from relativistic spin-orbit coupling 
effects, 5 ' 9,10 and treating finite-temperature effects, to name only two. 

Previous theoretical studies regarding the spin-wave dispersion in thin layers have been 
performed by Costa et al. 7,11 . The method used to investigate these systems consists of an 
empirical tight binding description, either from empirical data or from first principles calcu- 
lations, of both the magnetic layers and the substrate followed by a dynamical theory of the 
spin- waves. 12-14 Recently the formalism of atomistic spin-dynamics (ASD) was developed 15 
and was subsequently implemented in the UppASD package ' . With this technique it is 
possible to model spin waves at finite temperature, as shown recently for 1 ML of Fe on 
W(110). s The method relies on an accurate estimation of the exchange parameters in the 
Hamiltonian using first-principles calculations, and subsequently evolving it in time accord- 
ing to a modified Landau-Lifshitz-Gilbert equation of motion, 15 ' 18 as described in the next 
section. 

In this article, we extend the scope of our previous investigations, and employ the 
first-principles/ ASD strategy to calculate the magnon spectra of Co/Cu(001), Co/Cu(lll), 
Fe/Cu(001) and Fe/W(110), for varying thicknesses of the magnetic layer. The chosen sys- 
tems have been extensively studied experimentally and recently, using SPEELS, also magnon 
spectra have been measured for selected systems with certain thicknesses of the magnetic 
over layer. 
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I. THEORETICAL APPROACH 



First Principles Calculations and Atomistic Spin Dynamics Simulations 

Our strategy is to combine first-principles calculations with ASD simulations. In practice, 
we map an itinerant electron system onto an effective Hamiltonian with classical spins, 

W = "Hex + "Hso- 
This generalized Hamiltonian consists of an exchange term: 

Hex = ~ ^ Jjj^i • nxj, (1) 

where i and j are atomic indices, nij is the classical atomic moment and is the strength 
of the exchange interaction, and a spin-orbit interaction term: 

"Hso = K ^ ( m i ■ e K) 2 - ^2 Di J ( mi x m ^ ' ( 2 ) 

i i^j 

where K is the strength of the anisotropy field along the direction and D^- is the 
Dzyaloshiskii-Moriya (DM) vector. 19 ' 20 These represent the main relativistic effects present 
in low-dimensional magnets, although the magnitude of the DM term is often negligibly 
small. 

The ASD simulations rely on an accurate determination of the exchange integrals and 
spin-orbit coupling present in the generalized Hamiltonian. In order to determined these, 
we employed two different methodes for our first principles calculations. First, we set up a 
slab with 9-11 layers representing the substrate, having on one side the magnetic over- layer 
and a large region of vacuum, simulating the surface. The interlayer distances were relaxed 
according to the forces using the projector augmented wave method (PAW) as implemented 
in the VASP program 21 ' 2- '. Once the relaxed geometry is established, we employed the 
SPR-KKR package 23 ' 24 for calculations of the parameters to be used within the generalized 
Hamiltonian using the local spin density approximation (LSDA) for the exchange-correlation 
potential. The basis set consist of s, p, d and / orbitals and all relativistic effects are included 
by solving the Dirac equation. 

The exchange parameters Jjj and Dzyaloshinskii-Moriya (DM) vectors are obtained 
from the relativistic generalization by Udvardi et. a/ 6 ' 25 of the real-space method of in- 
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finitesimal rotations of Liechtenstein, Katsnelson and Gubanov (LKG), 26 ' 27 . Here, the full 
exchange tensor is calculated which can be decomposed into its isotropic part and its 
anti-symmetric part, which correponds to the DM- vector. The whole procedure relies on an 
adiabatic approximation in which the slow motion of the spins is decoupled from the fast 
motion of the itinerant electrons, a situation that is justified at low energy scales and for 
systems with reasonably large exchange splitting. 

Once calculated, the exchange parameters may be Fourier transformed to obtain the 
so-called adiabatic magnon spectrum. In the simple case of a single layer (corresponding to 
one atom per cell), the energy of a spin wave with respect to a ferromagnetic ground state 
is given by 



(3) 



where is the relative position vector connecting sites % and j. From this it is straight- 
forward to calculate the spin wave dispersion a;(q). 2s For systems with more than one atom 
per cell, as is the case for thin films consisting of more than one monolayer, the spin wave 
energies are given by the eigenvalues of the general N x N matrix here expressed in block 
form 



Ef^'-^(q) 



(4) 



••/"(<!)• Ef JS-J'K* 

where N is the number of atoms per cell (i.e. in this case the number of monolayers). 

Using in the ASD simulations the generalized Hamiltonian H as a starting point, the 
effective interaction field experienced by each atomic moment rrij is given by 

dU 



B, 



dnii 



The temporal evolution of the atomic spins at finite temperature is governed by Langevin 
dynamics, through coupled stochastic differential equations of the Landau-Lifshitz form, 



duii 



7 



■nii x [B; + bi(t)] - 7- 



a 



-m; x {nii x [Bi + b^t)]} 



(5) 



dt (l + a 2 ) 1 iy n 'm(l + a\ 

where 7 is the the electron gyromagnetic ratio. Temperature fluctuations are included via 
a stochastic Gaussian shaped magnetic field bj(i) with properties (bj(t)) = and 



{b k l (t)b l At')) = 2D8 l3 8 kl 5(t-t'), D 



a k B T 
(1 + a 2 ) iib^ 



(6) 



where i and j denote lattice sites, k and I the carteisian components and a is the Gilbert 
damping parameter which eventually brings the system to thermal equilibrium. 

The coupled equations of motion (5) can be viewed as describing the precession of each 
spin about an effective interaction field, with complexity arising from the fact that, since all 
spins are moving, the effective field is not static. In our calculations we evolve the stochastic 
Landau-Lifshitz equations using a semi-implicit method introduced by Mentink et a/. 18 

The principal advantage of combining first-principles calculations with the ASD approach 
is that it allows us to address the dynamical properties of spin systems at finite temper- 
atures. 16 ' 29 ' M> We focus in particular on two important quantities, the space- and time- 
displaced correlation function: 



where the angular brackets signify an ensemble average and k the cartesian component, and 
its Fourier Transform, the dynamical structure factor: 



where q and u are the momentum and energy transfer, respectively. The dynamical structure 
factor, S(q,u) is the quantity probed in neutron scattering experiments of bulk systems 31 , 
and can analogously be applied to SPEELS measurements. By plotting the peak positions of 
the structure factor along particular directions in reciprocal space, the spin wave dispersions 

may be obtained. 8 ' 16 ' 29 ' 30 

II. RESULTS 

A. Critical temperatures 

Within the Heisenberg model and in the case where the anisotropy and finite size effects 
are absent, a two dimensional magnetic film cannot sustain long range magnetic order at 




(7) 
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finite temperatures (the celebrated Mermin -Wagner theorem). However, the presence of 
symmetry breaking perturbations, such as small anisotropics, will overrule Mermin- Wagner 
theorem and finite temperature magnetism is possible, even in the case of a single magnetic 
layer. The critical temperatures are often reduced in thin films compared to their bulk 
counterpart. Using our ab initio calculated exchange interaction parameters, anisotropy 
and magnetic moments, we estimate the critical temperature, T c , in several thin magnetic 
films by means of Monte Carlo simulations and using the cumulant crossing method '". The 
systems we investigated are Cojv/Cu(111), Cojv/Cu(001), FeAr/Cu(001) and Fejv/W(110), 
where N indicates the number of magnetic atomic layers present. The results are summarized 
in Table I. It is difficult to directly compare the calculated T c values with experimental 
results 33,34 , due to the lack of information about interface quality, inter-diffusion between 
the substrate and the magnetic film and the surface roughness of the samples. All these 
imperfections affect the value of the critical temperature. 3 ' The calculations assume the 
idealized case with perfect interfaces and are likely to give the highest possible value of T c 
that the system can reach. Our data for 3 ML Co/Cu(lll) are in rather good agreement 
with the measured values for 2.5 ML Co/Cu(lll), although a comparison should be done 
with care, since the experiment and theory have slighly different Co thickness. For the 
1ML Co/Cu(lll) system as well as for the bulk Co, the theoretical values seem to be in 
agreement with measurements, whereas for the 2ML Co/Cu(001) system the calculated 
value of T c is two times larger than the measured value. For 2ML Fe/W(110) as well as 
for bulk Fe, theory and experiments are in good agreement, whereas for the Fe/Cu(001) 
system, theory consistently overestimates the measured values. Our calculations do however 
agree reasonably well with the results obtained by Pajda et. al using the random phase 
approximation (RPA). 
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B. Spin wave spectra 
1. Co on Cu(OOl) 

The fee Co/Cu(001) system represents a model system for the study of magnetic phe- 
nomena in thin films, since it does not exhibit strong structural, chemical or magnetic 
instabilities. 33 ' 42 It was therefore a natural candidate for investigation with SPEELS, and 
became the first ultrathin system in which its magnon spectrum was measured using this 
method. 2 The system could be described to a good level of accuracy within the context of 
the nearest-neighbour Heisenberg model on a semi-infinite substrate. In this case a surface 



Table I. Calculated values of the Curie temperatures Tq for Co and Fe overlayers on W and Cu 
substrates with different orientations and compared with published values determined theoretically 
and experimentally. 



System 


Theory: this work 


Tc (K) 

previous studies 


Experiment 


bulk fee Co 


1120 


1080 s7 , 125 38 , 1311 ± 4 39 


1388-1398 


Co N /Culll 








1 


255 




207-434 33 


2 


660 




500 [1.7 ML] 


Co N /Cu001 








1 


370 


426 36 




2 


680 




325 33 


3 


705 




650 [2.5 ML] 


8 


1070 






Fe N /Cu001 








1 


400 


515 36 


«100 


2 


380 




«270 


3 


560 




«290 33 


Fe N /W110 








2 


520 




450 41 
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mode exists with a dispersion curve hui^ sur ^ = 875(1 — cos(gao)) along the (110) direction, 
where J is the exchange coupling, S is the magnitude of the spin per primitive unit cell, 
q is the length of the magnon wave vector and ao = 2.55 A. The fit of this curve to the 
measured SPEELS data gave JS = 15.0 ±0.1 meV, which compares well with the value of 
JS = 14.7 ± 1.5 meV obtained from neutron data of bulk fee Co at long wave lengths 1 ' (i.e. 
in the regime q < 0.3 A -1 ). 

The Halle group followed up their initial report for the 8 ML film with results relevant 
to Co/Cu(001) systems with decreasing thickness down to 2.5 ML. 44 ' 45 These experiments 
demonstrated a very weak reduction in the energies required to excite the spin waves, relative 
to the bulk. 1. However, for all the thicknesses reported, the surface mode at the surface 
Brillouin zone boundary (the X point) is well below the bulk band edge, at around 240 meV. 
This difference is for the most part caused by the reduced number of nearest neighbours (NN) 
at the surface (only 8 atoms) with respect to the bulk case, where there are 12 NN present. 2 
For energy ranges where Stoner excitations are important, no direct comparison between 
calculated and experimental data for bulk Co is possible at the moment. 

Theoretical investigations of Co/Cu(001) thin films have been carried out using both the 
adiabatic approximation' 1 ' 36 and the random phase approximation (RPA) to a description 
of the spin response of the itinerant electron system. 12,13 The primary difference between 
the two approaches concerns the treatment of the particle hole excitations known as Stoner 
excitations, 15 which are neglected in the adiabatic case. Broadly speaking, these are not 
relevant at low energies, where the density of magnon states dominates. However, at higher 
energies or large wave- vectors their role becomes significant, and one would therefore expect 
theories based on the adiabatic approach to break down at the Brillouin zone boundary. In 
addition, Stoner excitations are more pronounced for low dimensional magnets compared to 
bulk magnets. 

Fig. 1(a) displays the calculated adiabatic magnon spectrum obtained for 8 ML Co/Cu(001), 
using Eq. (4) with values obtained using the LKG method. ' : The most notable fea- 
ture is the presence of several branches, one for each Co layer present. This is in contrast 
to experimental observations, where only the lowest ("acoustic") branch (filled circles in 
Fig 1(a)) as well as the second lowest branch 11 are observed. As Vollmer et al. point out, 2 
this indicates the shortcomings of a direct interpretation of their data in terms of a naive 
Heisenberg model. 
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Fig. 1(b) displays the adiabatic magnon spectrum obtained for 3 ML Fe/Cu(001). Un- 
fortunately, there are no experimental data to compare with. The values for the spin wave 
stiffness, D, can be estimated from the calculated adiabatic magnon spectra by measuring 
the curvature of the dispersions as q — » 0. For this system we find D being approximately 
210 meV A 2 . For the 1 ML Co on Cu(001) case we obtain a value of the order of 420 meV 
A 2 - a 15 % overestimate of the experimentally determined value of 360 meV A 2 , but con- 
siderably softer than the theoretical value determined by Pajda et al. of 532 ± 9, obtained 
by a real-space adiabatic approach. For 1 ML Fe on Cu(001) we obtain a lower value, of 
260 meV A 2 compared to the value of 331 meV A 2 by Pajda et a/. 36 . 




(a) (b) 

Figure 1. : Adiabatic magnon spectra (full lines) obtained from SPR-KKR calculations for (a) 8 
ML Co/Cu(001) together with experimental SPEELS data (circles) and (b) 3 ML Fe/Cu(001). 



As explained above, the adiabatic approximation becomes questionable at higher energies 
and/or large wave-vecors, since Stoner excitations become relevant. In order to address this 
issue, Costa et al. 12 ' 13 have developed a theory that explicitly takes these excitations into 
account. Their approach successfully describes the SPEELS measurements for the 8 ML 
Co/Cu(001) system: their calculated magnon spectrum is in agreement with the experiment 
over all the T — X line and correctly predicts a broadening of the "acoustic" spin wave peaks, 
along with an absence of standing spin waves giving rise to "optical" branches. Their work 
strongly indicates that a process known as Landau damping, through which the spin waves 
decay into the Stoner continuum, is at play in the Co/Cu(001) system. 

At present, the ASD method cannot handle one-particle Stoner excitations directly (these 
are responsible for the longitudinal fluctuations of the magnetic moments), although recently 
efforts have been made to include these effects in the LLG-equation, both directly 1 ' and 
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indirectly . However, the ASD method allows to address an alternative source of damping, 
namely the temperature, and we have performed simulations for different thicknesses of Co 
overlayers (1, 2, 3 and 8 atomic layers) under two very different conditions (Figs. 2 - 6). 
First we performed simulations at low temperature (1 K), where the damping constant a in 
Eq. (5) is set to 3 x 10 ~ 4 , a value that ensures a very weak coupling to the temperature 
bath. In this case the temperature effects are deliberately kept to a minimum. Secondly, we 
perform simulations in more realistic conditions, namely at room temperature (T=300 K) 
and a physically more plausible damping constant a = 0.05. All these results are discussed 
in detail in the following sections. 

1 ML 

Fig. 2 displays the spin wave spectra in a 1 ML thick Co layer on Cu(001). In the 
low temperature case (Fig. 2(a)), the spectra obtained from the atomistic spin dynamics 
simulations is almost indistinguishable from the adiabatic spectra. We point out that these 
results are obtained from a calculation of the dynamical structure factor S*(q, u). In order 
to facilitate the comparison between the spin wave spectra for a single layer of Co on Cu for 
different orientations of the substrate, we performed the calculations for both Cu(001) and 
Cu(lll) orientations (discussed in detailed below), at the same value of the temperature, 
i.e. T= 200 K. Increasing the temperature and damping (Fig. 2(b)) causes some broadening 
of the dynamical structure factor and as a result the magnon energy at the BZ boundary (X 
and M points) decreases slightly. However, overall the effects of temperature and dynamic 
treatment are not significant in this case. 

2 and 3 ML 

Qualitatively, the 2 ML and 3 ML case are not so different so we only display the cal- 
culated dynamic structure factor for 2 ML Co in Figure 3. Having more than one layer on 
top of Cu(001), we expect, apart from the acoustic branch, a number of optical branches 
to appear in the magnon spectra, similar to what is found in the adiabatic spectra. In the 
case of 2 ML we expect one branch of each kind, but as noticed in Fig. 3(a), even at very 
low temperature (T= 1 K) and extremely small damping, the optical branch is very weak, 
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(a) (b) 

Figure 2. Spin wave dispersion spectra obtained from ASD simulations of 1 ML Co/Cu(001) at 
(a) T= IK and small damping constant a = 3 x 10 _4 and (b) T = 200 K and realistic damping 
constant a = 0.05 

especially at small wave vectors, close to the T-point. In Section II C below, we make a 
deeper analysis of the dynamical structure factor to better understand the suppression of 
the optical branches. This type of analysis has been recently published by Taroni et al.. 49 
Using a more reasonable temperature and damping (Fig. 3(b)) erases most traces of the 
optical branch which is smeared out and remains slightly visible only in the X — M region. 
The magnon energy at the BZ boundaries is reduced in this case by roughly 25 meV which 
can be attributed to temperature effects. 

8 ML 

Fig. 4 displays the spin wave spectra on the full 8 ML Co stack on Cu(001). It is 
immediately noticed that the contrast between the calculated spectrum based on the dy- 
namical structure factor (Fig. 4(a)) with the spectrum directly obtained from first-principles 
(Fig. 1(a)) is strong: for values of q close to the T-point. In Fig. 4, there is little trace of the 
"optical" branches even in the case of low temperature and damping. We also note that at 
300 K, the magnon energies at the zone boundaries are reduced by roughly 25 meV, similar 
to the 2 ML case. The suppression of optical branches arising from a dynamical treatment 
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Figure 3. Spin wave dispersion spectra obtained from ASD simulations of 2 ML Co/Cu(001) at 
a) T= IK and small damping constant a = 3 x 10 ~ 4 and b) T = 300 K and realistic damping 
constant a = 0.05 

has also been pointed out by Costa et al. 12 ' 13 . We also note that the agreement between the 
theoretical acoustic branch that we calculated and experiment is rather good. 



2. Co on Cu(lll) 

For completeness, we also perfomed simulations on a single layer of Co on the Cu(lll) 
surface. Due to the fact that the critical temperature of the single Co layer on Cu(lll) is 
lower than the room temperature (i.e. T c =255 K), for this system we performed the high- 
temperature calculations at 200 K instead of room temperature. In Fig. 5, the spin wave 
spectra obtained at T=200K is displayed. The calculated spin wave spectra is noted to be 
softer than that of the 001-surface, which is consistent with the lower Curie temperatures 
obtained for this system. The calculated spin wave stiffness has a value of around 220 
meVA 2 , almost half of the value compare to what was found for the 001-surface. 
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(a) (b) 

Figure 4. Spin wave dispersion spectra obtained from ASD simulations of 8 ML Co/Cu(001) at 
a) T= IK and small damping constant a = 3 x 10 ~ 4 and b) T = 300 K and realistic damping 
constant a = 0.05. The experimental values are marked by white squares. 




Iql 



Figure 5. Spin wave dispersion spectra obtained from ASD simulations of 1 ML Co/Cu(lll) at 
T = 200 K and realistic damping constant a = 0.05 

3. Fe on Cu(OOl) 

The Fe/Cu(001) system is more complex than Co/Cu(001) for a number of reasons. Since 
Fe doesn't naturally exist in the fee phase it is difficult to grow thicker layers of Fe with good 

13 



quality. It is also well known, from theory, that bulk fee Fe exhibits a complicated magnetic 
phase diagram with many magnetic configurations with similar energies. The same thing 
is true for thick Fe layers on Cu. However, for thin Fe layers (1 to 3 ML), it is generally 
accepted that Fe adopts a ferromagnetic configuration. In the case of thicker layers there 
are several proposed magnetic configurations, for instance Sandratskii claimed that the 
magnetic structure takes the form J/|T for the three upper layers. We performed Monte Carlo 
simulations using exchange parameters starting from either ferromagnetic or the proposed 
magnetic structure of Sandratskii as reference state for the 3 ML case. Regardless of the 
starting configuration, we always obtained the ferromagnetic configuration as the ground 
state magnetic structure also for 3 ML Fe on CuOOl. However, the spin waves (Fig. 6) are 
soft which is reflected in the calculated adiabatic spin wave stiffness constant D for which 
we obtain a value of 170 meVA 2 , which is considerably lower than what we obtain for 1 
ML. We would like to point out that in general ASD gives a more realistic description of 
the magnon spectra than the adiabatic approximation (Fig. 1(b)). 

In Fig. 6 we show the results of spin wave spectra only for the 3 ML case. The spectra is 
rather different from Co/Cu(001) and the most striking difference is the overall softness of 
the magnons. However, we should keep in mind that the critical temperature for the 3 ML 
Fe/Cu(001) is much lower than for the analogous Co system. The optical branch for the Fe 
system is more pronounced at low temperatures but similar to the Co case it is suppressed 
close to the f-point for reciprocal vectors inside the Brillouin zone. For realistic conditions, 
i.e. at room temperature and for larger damping, the intensity of the optical branches is 
very weak and smeared traces of those branches remains visible mostly in the X — M region. 

4. Fe on W(110) 

In a previous study' s , we performed calculations of the spin wave spectra of a single Fe layer 
on top of W(110) and compared it with the experimentally determined spectra 5 and found a 
good agreement. In particular, we observed a softening of the magnons with respect to the 
bulk spin wave dispersion spectra. This magnon softening in thin layers is expected even 
when considering a simple NN Heisenberg model for the description of the dispersion curves 
in ferromagnets, where the exchange parameters are weaker due to a lower coordination at 
the surface. Here we focus on the magnons of a Fe bilayer on W(110). Experimentally, the 
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(a) (b) 

Figure 6. Spin wave dispersion spectra obtained from ASD simulations of 3 ML Fe/Cu(001) at 
a) T= IK and small damping constant a = 3 x 10 ~ 4 and b) T = 300 K and realistic damping 
constant a = 0.05 

magnon dispersion has been measured by SPEELS ' 1 where in addition the magnon lifetimes 
were determined ' 2 . In Fig. 7 we show our calculated magnon dispersion along symmetry 
lines in the first Brilloiun zone together with experimental data. Along the line T — H, we 
find a very good agreement between experiment and theory (Fig. 7(a)). If the momentum 
transfer is changed, so that one goes outside the first Brillouin zone, the different spin wave 
branches can be determined. This is analyzed in detail in Section II C. In Fig. 7(b) we show 
the so obtained magnon dispersion, where the optical modes can be clearly observed. 

Udvardi and Szunyogh predicted the asymmetry in the magnon spectrum arising from 
the Dzyaloshinskii-Moriya interaction. Later on it was experimentally detected by Zakeri 
et al. in a Fe bilayer on W(110). The spin wave asymmetry AE is defined as the dif- 
ference in the spin wave energy w(q) between q and — q, i.e AE = u;(q) — u{— q). If 
Dzyaloshinskii-Moriya interactions are absent, i.e. very weak effect from spin-orbit interac- 
tion, the asymmetry is zero for every wave vector q. In Fig. 8, we show our calculated spin 
wave asymmetry for the Fe bilayer on W(110) for wave vectors ranging from — H to H in the 
two dimensional Brillouin zone, using theoretically determined Dzyaloshinskii-Moriya inter- 
action parameters. The simulations were performed at room temperature, as in experiment, 
with realistic damping. We obtain a qualitatively good agreement with experiment but 
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Figure 7. Spin wave dispersion spectra obtained from ASD simulations of 2 ML Fe/W(110) at 
T=300 K and a = 0.01. a) Sampling inside the first Brilloiun zone (r = [00]) b) Sampling shifted 
by vector r = [10]. Experimental values are marked by white squares. . See the text for details. 

the amplitude of the calculated spin wave asymmetry is slightly overestimated ( ~ YlmeV 
compared to ~ 8meV in experiment). There are several explainations to this discrepency, 
primarily the asymmetry is sensitive to the value of the Dzyaloshinskii-Moriya interaction 
which is very delicate to calculate from ab-initio theory. In particular, it requires many 
k-points in order to converge as well depends a bit on technical parameters such as basis-set, 
sphere radii and choice of exchange correlation functional. The assumption made in the cal- 
culations of an atomicaly sharp interface between Fe and W may also be a limiting factor. 

C. Analysis of the dynamical structure factor 

In neutron scattering, the susceptibility can be written in the following form ' 



where xa c and xop are the susceptibilites originating from accoustic and optical branch, 
respectively, q is a reciprocal vector within the primitive Brillouin zone (BZ), r = [hkl] = 
/i&i + kb 2 + lb 3 , (j> — t • p, where p is a vector connecting two sublattices. In this manner, 




(9) 
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Figure 8. (Color online) Calculated spin wave asymmetry for the magnon spectrum of 2 ML 
Fe/W(110), using theoretically determined Dzyaloshinskii-Moriya interactions. The experimental 
values have been obtained by Zakeri et al. for M || [110]. 

by changing the momentum transfer by varying r, the intensity of the accoustic and optical 
branches is changing. If we take the Fe bilayer on W(110) as an example, r = [0 ^ 0.5]a and 
the reciprocal vectors are restricted in the film-plane (/=0), then it follows that <p = fm + kn. 
Inside the primitive BZ, the phase = and in the limit q — > 0, in Eq. (9) the acoustic 
term will dominate and will be detected in experiment. If we go outside the first BZ, it is 
possible to have a situation where the optical term dominates, on the expense of the acoustic 
response, for instance by choosing r = [10], as illustrated in Fig. 7(b). If there are more 
than two atoms in the unit cell, the analysis becomes more complicated but the principle is 
the same. 

III. CONCLUSIONS 

Summarizing, we have shown that combining first principles calculations with atomistic 
spin dynamics simulations provides a powerful tool for studies of magnetic excitations in 
low dimensional systems. Our method can be applied to other low-dimensional systems, 
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particularly for cases in which exotic non-collinear ground states occur. 5 ' 9,10 It is hoped 
that the remarkable progress using the SPEELS technique will continue to provide new 
and surprising experimental results, for which the currently presented theory seems to be 
a good tool for analyzing the experimental dat and predicting magnon spectra and related 
properties. The materials presented here have several common features, e.g. the absence of 
optical modes in the magnon curves, both as determined by experiments as well as obtained 
by theory. This fact has been analyzed in detail and it is argued that also the optical modes 
should be visible if one considers excitations which allow for momentum transfer outside the 
first Brillouin zone. The realization of this in the SPEELS method is clearly a challenge. 
In addition, we show that all thin film systems investigated here have a spin wave stiffness 
which is considerably softer compared to the bulk value. This applies both to Fe as well 
as Co films. Finally, we report on a quantitative agreement between theory and measured 
magnon curves for 2 ML Fe on W(110) and 8 ML Co on Cu(OOl). 
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